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ABSTRACT 

We present observational constraints on the nature of dark energy using the Super- 



nova Legacy Survey three year sample (SNLS3) of Guy et al. (2010) and Conley et al 
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(2011). We use the 472 SNe la in this sample, accounting for recently discovered cor- 
relations between SN la luminosity and host galaxy properties, and include the effects 
of all identified systematic uncertainties directly in the cosmological fits. Combining 
the SNLS3 data with the full WMAP7 power spectrum, the Sloan Digital Sky Survey 
luminous red galaxy power spectrum, and a prior on the Hubble constant Hq from 
SHOES, in a flat universe we find flm = 0.269 ±0.015 and w = -1.06lt[J;[J^^ (where the 
uncertainties include all statistical and SN la systematic errors) - a 6.5% measure of 
the dark energy equation-of-state parameter w. The statistical and systematic uncer- 
tainties are approximately equal, with the systematic uncertainties dominated by the 
photometric calibration of the SN la fluxes - without these calibration effects, system- 
atics contribute only a ~ 2% error in w. When relaxing the assumption of flatness, 
we find Qm = 0.271 ± 0.015, Qk = -0.002 ± 0.006, and w = -1.069t[j:j]^^. Parameter- 
izing the time evolution of w as w{a) = wq + Wa{i — a), gives wq = —0.905 it 0.196, 
Wa = — 0.984]txo97 in a flat universe. All of our results are consistent with a flat, w = —1 
universe. 

The size of the SNLS3 sample allows various tests to be performed with the SNe 
segregated according to their light curve and host galaxy properties. We find that the 
cosmological constraints derived from these different sub-samples are consistent. There 
is evidence that the coefficient, /3, relating SN la luminosity and color, varies with host 
parameters at > 4(T significance (in addition to the known SN luminosity-host relation) ; 
however this has only a small effect on the cosmological results and is currently a sub- 
dominant systematic. 

Subject headings: cosmology: observations - cosmological parameters - supernovae: 
general - surveys 



1. Introduction 



The discovery of the accelerating universe ranks as one of science's landmark achievements 
in the 20th century. Surveys exploiting distant Type la supernovae (SNe la) as standardizable 



candles (Riess et al. 1998 Perlmutter et al. 1999) revealed the presence of a "dark energy" that 



opposes gravity and accelerates the expansion of the Universe. When these SN observations are 



combined with measures of large-scale structure (Cole et al. 2005; Eisenstein et al. 2005 Percival 



et al.||2007t |Reid et al.||2010D and the cosmic microwave background (CMB; e.g., |de Bernardis et al. 
2002t [Bennett et al.|2003t [Larson et al.|201lD , this dark energy emerges as the dominant component 
of the Universe responsible for 70-75% of its energy density at the present epoch. 

A compelling physical explanation of dark energy remains distant despite a range of possibil- 



ities being postulated (for reviews see Copeland et al. 2006 Frieman et al. 2008). Astrophysical 



measurements of the dark energy's equation-of-state parameter w (the ratio of its pressure to den- 
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sity, p/p), and its variation over cosmic history, can help distinguish the possibihties. The classical 
"Cosmological Constant" is equivalent to a vacuum energy density with negative pressure, constant 
in time and space: w = —1. The broad family of "quintessence" models, a dynamical form of scalar 
energy field, mostly predict —l<w<—^. A measurement of w < —1 would be a signal of even 
more exotic physics. 



SNe la remain, at present, the most direct and mature method of probing this dark energy due 
to several decades of intensive study and use in cosmology (see the review of |Howell|20T0| ). Thought 
to be the result of the thermonuclear destruction of an accreting CO white dwarf star approaching 
the Chandrasekhar mass limit (e.g., Hillebrandt Sz Niemeyer|2000 ) , they are standardizable candles 
which explode with nearly the same brightness everywhere in the Universe due to the uniformity 
of the triggering mass and hence the available nuclear fuel. Their cosmological use exploits simple 
empirical relations between their luminosity and other parameters. Brighter SNe la have slower. 



wider light curves (Phillips 1993) and are typically bluer than their faster, fainter counterparts 
(Riess et al. 1996 Tripp 1998). Detailed searches for a "third parameter" have recently shown 
that, after applying these first two corrections, brighter SNe la reside in more massive host galaxies 
dKelly et al.]|2010[ [Sullivan et ar]|2010[ [Lampeitl et"ar]|2010D . 

The application of relations between SN la light-curve-shape, color, and host galaxy proper- 
ties provides robust distance estimates which allow SNe la to be used to measure cosmological 



parameters. When combined together, recent SN la samples ( 


Astier et al. 


2006 


Riess et al. 


2007 


Miknaitis et al. 


2007; 


Holtzman et al. 


2008 


Hicken et al. 


2009a 


Contreras et al. 


2010) provide 



measures of dark energy generally consistent with a cosmological constant of w = —1 with statis- 
tical uncertainties in w of around 5-7%, and systematic uncertainties of 8-14%, depending on the 



method used and assumptions made (Kessler et al. 2009 Hicken et al. 2009b Amanullah et al. 



2010). These SN la samples are now sufficiently large that our understanding of systematic uncer- 



tainties has a direct impact on our measurement of dark energy (Kessler et al. 2009: Conley et al. 



2011), particularly when combining SNe la from many different surveys at different observatories. 



Systematic uncertainties which affect the cosmological analysis of SNe la arise from two broad 
sources. The first is experimental systematics; for example, photometric calibration or survey selec- 
tion biases. Due to the correlations between SN la magnitudes that these uncertainties introduce, 
accounting for their effects in the cosmological fits is particularly important. Fortunately this is a 
tractable problem as the uncertainties are understood in modern SN la samples and can be accu- 



rately modeled, albeit only after detailed and painstaking work (Regnault et al. 2009 Guy et al. 
20T0l [Conley et al.|[20lT| ). The second source of systematics arises from an incomplete understand- 
ing of their astrophysics (e.g., progenitor configuration, mass transfer and explosion mechanism, 
etc.). The most pernicious possibilities include evolution in SN la properties with redshift tracking 
changing metallicities/ages of the progenitor stars, and varying dust extinction or color laws; the 



correct treatment of SN la color-luminosity relationships are particularly uncertain (e.g., Conley 



et al. 2007). The effects of these potential systematics are more nebulous due to the difficulty in 



modeling SNe la explosions, but can be investigated empirically. Studies which compare local SN la 
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spectra with those at high-redshift find a remarkable degree of similarity across 2780A to 6OOOA 
( Hook et al.]|2005t|Blondin et al.|2006l [Foley et al.|2008tpiis et al.||2008HBalland et al.|2009t[Cooke 
et al.|[20lT ) with only small differences in the strengths of some intermediate mass element features 



(Sullivan et al. 2009), consistent with a mildly evolving mix in SN la demographics with redshift 



expected from popular SN la delay-time distribution models (Mannucci et al.||2005 Sullivan et al 



2006b; Mannucci et al.||2006 Howell et al.||2007 ). To date, no definitive evolutionary signature with 
redshift, which would directly impact a cosmological analysis, has been located. 

This paper presents the cosmological analysis of the three year Supernova Legacy Survey SN la 
sample (SNLS3). Our sample and methods are presented over the course of several papers. The 



first, Regnault et al. (2009), deals with the photometric calibration of the SNLS SN la fluxes 



and associated systematic uncertainties, including corrections for spatial non-uniformities in the 



SNLS photometric imager. Guy et al. (2010, hereafter GIG), presents the light curves of the SNLS 
SNe la themselves, together with a comparison of SN light curve fitting techniques, color laws. 



systematics and parameterizations. Conley et al. (2011, hereafter Gil), discusses systematic effects 



in the cosmological analysis (including covariance matrices accounting for correlations between the 
distances to different SNe), presents light curve parameterizations of external SNe la used in the 
analysis, describes the various light curve quality cuts made to produce the combined sample, and 
provides the cosmological constraints obtained from the SN la data alone. Other papers describe 



the SNLS selection biases (Perrett et al. 2010), the SN la host galaxy information (Sullivan et al. 



2010), and the spectroscopic confirmation and redshift measurements ( Howell et al.|[2005 Bronder 



et al. 2008 Balland et al. 2009: Walker et al. 2011). This paper performs a cosmological analysis 



combining the SN-only analysis of Gil with other external, non-SNe constraints. 



We use 242 well-sampled SNe la over 0.08 < z < 1.06 from the SNLS together with a large 
literature sample: 123 SNe la at low-redshift, 14 SNe la at 2; > 0.8 from the Hubble Space Telescope, 
and 93 SNe la at intermediate redshift from the first year of the Sloan Digital Sky Survey-H SN 
search. We include the effects of identified systematic uncertainties directly in our cosmological 



fitting analysis using an approach outlined in detail in Gil This allows our cosmological parameter 



uncertainties to include systematic as well as statistical uncertainties, with covariances between 
different SNe which infiuence the cosmological fits accounted for. Examples of effects which cause 
such covariances include common photometric zeropoints for different SNe, or selection effects for 
SNe from the same survey. Appropriate covariance matrices allowing other users of this combined 



data set to directly include systematic effects in subsequent analyses can be found in Gil 



The advantages of the enlarged SNLS data set are multiple. Most obviously, this represents a 
threefold increase in the SNLS sample size compared to the first year SNLS cosmological analysis 



presented in Astier et al. (2006 hereafter A06), and as such provides a significant improvement in 



the statistical precision of the cosmological constraints. Several improvements in survey strategy 
were made following the first year of SNLS, including a more regular observing cadence together 
with longer z-band exposures, important for the highest redshift events. Moreover, the enlarged 
data set allows sources of potential astrophysical systematics to be examined by dividing our SN la 
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sample according to properties of either the SN (e.g., hght curve width) or its environment (Sulhvan 



et al. 2010). The increased size of the SNLS data set has also enabled a better understanding 



of SN la light curve and spectral properties (particularly at A < 3600 A in the rest-frame) with a 



corresponding improvement in the methods for estimating their distances (Hsiao et al. 1 20071 Guy 



et al. 


2007 


Ellis et al. 


2008 


Conley et al. 


2008 



years of the SNLS data also allow an improved photometric calibration of the light curves and a 
more consistent understanding of the experimental characteristics ( Regnault et al.||2009| ). 



A plan of the paper follows. §[2]provides a brief overview of the SN la data, and §[3]describes our 
methodology for determining the cosmological parameters. Our cosmological results are presented 
in § |4j § [5] discusses cosmological fits to various sub-samples of our SN population designed to 
assess possibilities of astrophysical biases within the SN la sample. We summarize and conclude in 

m 



2. Supernova Data and Methodology Overview 



We begin by briefly reviewing the SN la datasets and the various techniques that we use in 
the cosmological analysis. Full details of all of our procedures can be found in |G10 and Cll 



as 



well as Guy et al. (2007), Conley et al. (2008), Perrett et al. (2010), and Sullivan et al. (2010). 



2.1. The SN la samples 



Our SN la samples are divided into two categories: those discovered and confirmed by the 
Supernova Legacy Survey (SNLS), and those taken from the literature which sample different 
redshift ranges to SNLS. The SNLS uses data taken as part of the five-year Canada- France-Hawaii 
Telescope Legacy Survey (CFHT-LS). CFHT-LS is an optical imaging survey, the deep component 
of which conducted repeat imaging of 4 fields every 3-4 nights in dark time with four filters, allowing 



the construction of high-quality multi-color SN light curves (GIO). Spectroscopic follow-up is used 



to confirm SN types and measure redshifts, critical in obtaining clean samples of SNe la as reliable 



photometric identification techniques have yet to be developed, despite recent progress (Kessler 



et al. 2010 Bazin et al. 2011). Candidates were prioritized following the procedure outlined in 
Sullivan et al. (2006a). SNLS benefited from large time allocations on 8-lOm class telescopes - 
~ 1500 hours over five years - including the Gemini North and South telescopes, the European 
Southern Observatory Very Large Telescopes, and the Keck telescopes. Nearly all our SN spectra 



Bronder et al.|2008 


Elhs et al.|2008 


Balland et al.|2009 


Walker 



et al.|[20lT| >. All spectra are analyzed and uniformly typed according to the classification schemes 



^Spectra and light curves for the SNLS3 sample are available at the University of Toronto's Research Repository, 
T-Space: jhttps : //tspace . libraLry .utoronto . ca/snl s, as well as in the cited papers. 
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oflHowell et al.l (|2005|) andlBalland et al.|(|2009p. Further information on all 242 SNLS SNe la that 



we use, including light curve parameterizations, can be found in GIO 



The SNLS dataset is complemented with SNe la from the literature over redshift ranges that 
the SNLS sample does not cover. We use 123 SNe la at low-redshift {z < 0.08) from a variety 



of sources (primarily Hamuy et al. 1996 Riess et al. 1999 Jha et al. 2006 Hicken et al. 2009a 



Hicken et al. 


2009a 


of 


Riess et al. 


(2007), 



and 93 SNe la over 0.06 ^ z < 0.4 from the first year of the Sloan Digital Sky Survey-II SN search 
( Holtzman et al.|2008 ). Light curve parameterizations and other data for these events, on the same 



relative system as that of GIO, can be found in Cll 



We also considered including other SN la samples that, at least in part, probe the same redshift 



range as SNLS. However, we do not do this for several reasons (see also Cll for a discussion of 
these points). SNLS is designed to control systematics as much as is possible - a single telescope 
survey with a well-understood photometric response and calibration, using deep exposures in fil- 
ters that allow the same rest-frame colors to be measured for most of the redshift range that it 



probes. There is also significant published information on the host galaxies (Sullivan et al. 2010), 
essential for the cosmological analysis and which is not available for other higher-redshift samples. 
Finally, SNLS is by far the largest and best observed (i.e., highest signal-to-noise for each event) 
SN la sample over 0.3 < z < 1.0. Adding other SNe to this might lead to marginally improved 
cosmological constraints from a purely statistical perspective, but would certainly lead to a much 
more complex and uncertain analysis of systematic uncertainties when combining data from many 
surveys conducted at many telescopes. 



All SNe are corrected for Galactic extinction, Malmquist and other selection biases (Perrett 



et al. 2010; Conley et al. 2011), and peculiar velocities (at low redshift). 



2.2. Light curve fitting and distance estimation 



We parameterize the SN la light curves for distance estimation using the combined results 
from updated versions of two independent light curve fitters - SiFTO (Conley et al. 2008) and 



SALT2 (Guy et al. 2007). Both techniques provide an estimate of the SN peak rest-frame i?-band 



apparent magnitude at the epoch of maximum light in that filter, a measure of the SN light curve 
shape, and an estimate of the SN optical B — V color (C). SiFTO parameterizes the light curve 
in terms of stretch (s), while SALT2 uses a related parameter xi. The two light curve fitters are 
compared in G10[ which also provides details of the techniques used to average them into a single 
light curve parameterization for subsequent distance estimation. The distance estimation technique 
used, based on the combined SiFT0/SALT2 light curve parameters, is described in § [sj 

Other light curve fitting and distance estimation techniques are available. In particular, the 
MLCS2k2 fitter ( |Jha et al.|2007l ) has been widely used in previous SN la analyses (e.g.lHicken et al. 
2009bt [Kessler et ar]|2009[ ), and the use of MLCS2k2 versus SALT2 has led to significantly different 
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cosmological parameters in some cases. The two techniques are mathematically equivalent (to first 
order - see section 4.2.3 in GIO), and many of the apparent differences can instead be traced back 



to the training data and priors on SN color. We do not use MLCS2k2 in the SNLS3 analysis for 



several practical reasons, discussed at length in GIO and Cll, and which we summarize here 



There are apparent calibration problems with observer [/-band SN la data (which we do not 
use anywhere in our analysis) which MLCS2k2 is reliant upon for its distance-estimation training. 
These include: 



• Observer-frame [/-band SN la data show more scatter around individual SN la light curve fits 
than can be accounted for by the published observational uncertainites. This large intrinsic 
scatter is not seen in SNLS and SDSS-SN SN observations transformed into rest-frame U- 
band. 

• In U — B versus B — V color-color space, the SNe with observer-frame [/-band data show a 
systematic offset compared with the other SN samples (see Fig. [T]), 

• There is significant tension between SNe with observer- frame [/-band data, and those without, 
in the cosmological fits. This tension disappears if this [/-band data is removed from the light 
curve fits for the low-z SNe. 



A sample of low-z SNe from the Carnegie Supernova Project ( Contreras et al.pOlO ) with observer- 
frame u' data show none of these three problems above. That is, the properties of the low-z SN la 
?7-band data are inconsistent with SNLS at z ~ 0.5, with SDSS-SN at z ~ 0.25, and with the CSP 
sample at z ~ 0.02. 



Various possibilities for the origin of this J7-band anomaly are outlined in Cll section 2.6 (see 
also Kessler et al.|[2009 ). For this to be an evolutionary or astrophysical effect, any evolution must 
be astonishingly sudden (i.e., turning on at z ~ 0.25 and then not evolving any further out to 
z = 1) - yet somehow only effect one of the two sources of low-z data (those with U, instead of u' , 
data). It must also somehow not manifest itself in maximum-light spectral comparisons between 



Foley et al. 


2008 


Cooke et al. 


2011 


is seen in Cooke et al. 


(2011 


is in 



that implied by the [/-band photometry problem. 

The overwhelming likelihood is that this is a problem with the notoriously difficult calibration 
of the observer- frame [/-band data - there is no evidence that it is an astrophysical effect, and 
significant evidence that it is not. Note that an MLCS2k2 trained without the [/-band data severely 
impacts the science that can be done - as a distance estimator, MLCS2k2 requires z < 0.06 SNe 
la for the training, and therefore cannot be supplemented with SNLS/SDSS SN data sampling the 
rest-frame [/-band as is the case with SALT2 and SiFTO. 



Even with a version of MLCS2k2 not trained using the current [/-band data, additional prob- 
lems remain. These are discussed at length in GIO in their section 4.2. In that analysis, the authors 
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noted that MLSC2k2 requires the use of priors that color variation in SNe la is caused by dust 
extinction, which can introduce additional biases into the estimated distances - there is no evidence 
that color variation in SNe la is caused purely by dust extinction, and significant evidence that 



intrinsic SN properties make the story more complex (e.g., ( 


I!onley et al. 


2007 


] 


^asen et al. 


2009 


Foley & Kasen 


2011 


Chotard et al.||2011 


Nordin et al.||2011 


Maeda et al.||2011t 


)). Additionally, a 



Cardelli et al. ( 1989 )-like color-variation law for the SNe is assumed, which again is not supported 



by the data (Guy et al. 2007, 2010). 



Disentangling the effect of intrinsic color variation from color variation due to dust would 
require additional parameters in MLCS2k2, SALT2 or SiFTO which currently do not exist - at 
present no light curve fitter correctly disentangles the two effects. In particular, a weakness of the 
approach used in this paper (see § [s]) is to conflate intrinsic and dust effects into a single parameter, 
/3, during the distance estimation. However, the tests that we are able to perform - for example 
examining the evolution in this P parameter with redshift (figure 14 in CIO) - show no evidence 
for any significant systematic effects in the SALT2/SiFT0 distance estimation method. The case 
of maximal evolution in (3 supported by our data are included in our error budget (see 



2.5). 



In conclusion, the differences reported in Kessler et al. ( 2009 ) are now understood and can be 



traced back to a combination of the [/-band calibration and color priors - and should not, therefore, 
be considered as systematics. These above issues were a significant motivation for us to develop 
two independent light curve fitters within SNLS that differ conceptually in the way that color is 



handled ( Guy et al. 2007 Conley et al. 2008 Guy et al. 2010 ). Diff'erences in the light curves fits 
from these two codes are carried through as an uncertainty in our analysis - so our final quoted 
errors on the cosmological parameters fully include this effect. 



2.3. SN Selection 



We apply various selection cuts to the SN la samples designed to ensure an adequate wavelength 



and phase coverage in the light curve fits. These are described in detail in GIO and Oil, and 



essentially ensure that an accurate light curve width, rest-frame color, and peak rest-frame B- 



magnitude can be measured. Additionally, each SN must be spectroscopically confirmed (see Howell 



et al. 2005, for a discussion of our spectroscopic classification criteria), have a minimum CMB- 
frame redshift (zcmb) of 0.010 (after peculiar velocity correction), be spectroscopically normal, 
have a Galactic extinction of E^'Zy < 0.2, and be of normal stretch (0.7 < s < 1.3) and color 



(—0.25 < C < 0.25). We also remove 6 outliers on the Hubble diagram - see Cll for details 



2.4. Host galaxy parameterizations 

Recent analyses have found correlations between SN la luminosity and their host galaxies, even 
after the well-known fainter-faster and fainter-redder corrections have been made ( Kelly et al.|201() 
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Sullivan et al.|2010 Lampeitl et al.|2010 ). At the present time, it is not clear from an astrophysical 
perspective which host galaxy parameter is the optimal choice to use (as the third parameter) in 



the analysis. The observed effects could be due to progenitor metallicity (e.g., Kasen et al. 2009), 
progenitor age (e.g., Krueger et al.|[2010 ), or possibly some other parameter. Empirically, trends in 
SN la luminosity are seen with a variety of derived host galaxy parameters, including stellar mass 
(Mgteiiar); Star formation rate (SFR), and inferred stellar age. However many of these parameters 
are strongly correlated when derived from broad-band photometry available for the host galaxies. 

We use the host galaxy stellar mass as the third variable in our cosmological fitting. This has 
the advantage of being easiest to determine from sometimes limited host data, and shows significant 
trends with the SN la luminosities. We derive the host information by following the method in 



Sullivan et al. (2010), fitting the broad-band spectral energy distribution (SED) of the host galaxies 
using simple galaxy population synthesis models ( Le Borgne &l Rocca-Volmerange||2002 ) . The host 



galaxy information can be found in GIO and Cll for all the SNe used in our analysis. 



2.5. Systematic uncertainties 



We consider a variety of systematic uncertainties in our analysis, discussed in detail in Cll 



Details of the construction of the covariance matrices that encode this information can be found in 
that paper. For each systematic we estimate its size and adjust that variable in the light curve fits. 
All the light curves are refit, including a re-training of the light curve model where required, and 
the derived light curve parameters {uib, s, C) compared for each SN with and without the inclusion 



of the systematic. These differences are converted into a covariance matrix following Cll 



3. Cosmological fit methodology 

Having summarized the main features of our dataset, we now turn to the cosmological analysis. 
We write the a-s 

^,^^ (,„^-„ ) 

SNe Stat ' "int 

where we have omitted the covariance error matrix for clarity. Cstat is the identified statistical error 
and includes uncertainties in both itib and m^°'^, a-mt parameterizes the intrinsic dispersion of each 
SN sample (see below), and the sum is over the SNe la entering the fit. thb are the maximum- light 
SN rest-frame i?-band apparent magnitudes and m^"'^ are the model -B-band magnitudes for each 
SN given by 



mB° = 5 logio '^L {Zheh Zcmh, W, Qrn, ^DE, ^k) 

-a{s-l) + (3C + MB 



(2) 
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where w is the equation of state parameter of dark energy, firn and Ode sue the fractional energy 
densities of matter and dark energy (for w = —1, J^de = ^a); is the curvature constant, and a 
and /3 parameterize the s and C-luminosity relationships. Any linear variation between SN intrinsic 
color and s will be absorbed into the a term, z^ei is the heliocentric redshift used in the light curve 
fits. Vl is the c/Hq reduced luminosity distance with the c/Hq factor absorbed into Mb (here c 
is the speed of light and Hq the Hubble constant). Explicitly, Mb = Mb + 51ogio(c///o) + 25, 
where Mb is the rest-frame absolute magnitude of a SN la in the i?-band. Neither Ho nor Mb are 
assumed during the fitting process. 

We allow Mb vary as a function of host galaxy stellar mass (Mgteiiar) to account for relations 
between SN la brightness and host properties that are not corrected for via the standard s and 



C-luminosity relations following Sullivan et al. (2010). Explicitly, we fit for M^ in galaxies with 
Msteiiar < 10^° Mq, and M% when M^teiiar > 10^° Mq. We could, of course, allow the other nuisance 
parameters a and f3 to vary according to host type - we discuss this further in § [5} 

The statistical errors affecting each SN include the statistical error in from the light curve 
fit, the statistical error in m^°'^ (essentially aag and /3ac), a peculiar velocity error of 150 kms~^ 
after correction for a local bulk flow model, the error in Zhei projected into magnitude space, a 



10% uncertainty from Milky Way extinction corrections (Schlegel et al. 1998), a random scatter 



due to gravitational lensing following Jonsson et al. (2010) of fiicns = 0.0552;, and the covariances 



between s, C and niB for an individual SN (these parameters are correlated as they are determined 
from the same light curve data), crgtat is updated during the fits as a and (3 are altered. The 
(Tint term parameterizes the extra dispersion in me required to give a P^r degree of freedom 



(DOF) of one in the cosmological fits (e.g., Perlmutter et al. 1999). This "intrinsic" dispersion 



arises from unidentified sources of error in our analysis, as well as the imperfect nature of SNe la 
as standard candles. cJint may also include contributions from unidentified experimental errors and 
survey selection effects, and there is no a priori reason for dint to be the same from SN sample to 
SN sample; we allow a different ai^t for each sample, the values for which can be found in table 4 
of Cll These values are not varied in the fits, but the values are fixed to give a x^/DOF = 1 for 



each sample in the SN only cosmological fits of Conley et al. (2011). Note that more sophisticated 



statistical techniques for treating Uint and its uncertainty have been proposed (March et al. 2011). 



To include systematic errors we generalize eqn. ([T]) by constructing a covariance matrix C 
to replace the a terms. C is the combination of a systematics covariance matrix Cgyst and two 
covariance matrices containing statistical uncertainties: Cgtat which contains statistical errors from 
the SN model used in the light curve fit and which are therefore correlated between SNe, and 
Dstat) a purely diagonal covariance matrix generated from the statistical errors described above. 
We include both Cstat and Dgtat when performing fits based only on statistical errors. 



We then minimize the according to 



X 



n 

N 



niB 



m- 



mod 



1-1 



me 



m 



mod I 



(3) 
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This methodology allows the quoted uncertainties on the fit parameters to directly include system- 
atic errors, as well as correctly accounting for systematic and statistical uncertainties which induce 
correlations between different SNe and thus alter the position of the best-fit cosmological model. 



3.1. Fitting techniques 

We use three approache^to perform our cosmological fits. For relatively simple cosmological 
fits involving a small number of parameters, we use a grid technique that computes the of eqn. Q 
at every point converting into a probability via P oc exp (— 2^^), with the proportionality set by 
normalizing over the grid. The "nuisance parameters" a, /3 and Mb are marginalized over when 
generating confidence contours in the parameters of interest, and we report the expectation value of 
the marginalized parameters. Due to the (relatively) fast run-time, and the contour visualization, 
this fitting technique is particularly well-suited to analyzing the magnitude of the individual sources 
of systematic uncertainty in our analysis, which would be impractical with more complex and slower 
fitting approaches. 

The second approach is a minimization routine which simply reports the best-fit. The 
results of this technique should be close to the reported values from the grid marginalization, but 
should not be expected to agree exactly, and we provide both. Note that the fiint calculated by jCll 



is performed for the marginalization approach fits - when these a^t are used in the x^ minimization 
fits, a x^/DOF < 1 should be expected as the best-fit parameters from the marginalization fits will 
not lie at a minimum in x^- 



The third approach is the CosmoMC program ( [Lewis fcBridle|2002D , which uses a Markov- Chain 



Monte Carlo technique to explore cosmological parameter space. We use this approach for our main 
cosmological results. We made the following modifications to the May 2010 version of CosmoMC 
package to handle SNLS3 data: first, we properly marginalize over the SN nuisance parameters a 
and (3 rather than holding them fixed; second, we keep track of the difference between heliocentric 
and CMB frame redshifts, which enter into the luminosity distance differently, important for some 
of the lowest-^ SNe; third, we have added the ability to fit for the host-dependence of SN la absolute 



luminosities as described in Sullivan et al. (2010). 



The first item above is handled most efficiently by explicitly fitting for a and /3 along with 
the cosmological parameters, as internally marginalizing over their values is computationally more 
expensive for the SNLS3 sample. The consequences of incorrectly holding the nuisance variables 
fixed, or of simply substituting the values that minimize the x^i both true of the default CosmoMC 
SN la implementation, are discussed in section 4.6 and Appendix C of Cll, as well as § [4] of this 



■^All the computer programs and code referred to in this paper are available at https : //tspace . library . 
"utoronto . ca/snls along with the SN la light curves, spectra, light-curve parameters, covariance matrices, and 
some of the CosmoMC chains. 
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paper. In brief, this simplified approach leads to both underestimated uncertainties, and biased 
parameter estimates, due to small correlations between a, /3, and the cosmological parameters. 

For CosmoMC fits where we allow for a time- varying dark energy equation of state {w{a) = 



Wo + Wa{l — a), where a is the scale factor), we follow the prescription of Fang et al. (2008). 
Further, in the CosmoMC fits we do not consider massive neutrinos, and assume a simple power-law 
primordial power spectrum (i.e., we neglect tensor modes, and any running of the scalar spectral 
index) . 



3.2. External, non-SN datasets 



We include several external non-SN datasets in our fits. For the grid marginalization and 
minimization fits, we use two external constraints. The first is the Sloan Digital Sky Survey (SDSS) 



Data Release 7 (DR7) Baryon Acoustic Oscillations (BAO) measurements of Percival et al. (2010). 
This is a Gaussian prior on the distance ratios rs{zd)/ Dv{z) at z = 0.2 and z = 0.35, where rs{zii) is 
the comoving sound horizon at the baryon drag epoch, and Dy{z) is a spherically averaged effective 
distance measure given by Dv{z) = [{l + zfD\{z)cz/H{z)Yl^ ( [Eisenstein et al.|2005[ ) where Da{z) 
is the proper angular diameter distance. The second is a prior based on the Wilkinson Microwave 



Anisotropy Probe 7- year (WMAP7) "shift" parameter R (Bond et al. 1997), the "acoustic scale" 



la, and the decoupling redshift z*, as defined in Komatsu et al. (2011), following the prescription 



of Komatsu et al. ( 2009 ) . This prior includes most of the power of the CMB data for measuring 



dark energy (e.g. Wang et al. 2007). 



For our main cosmological fits, with the CosmoMC program, we use different external constraints: 
the power spectrum of luminous red galaxies (LRGs) in the SDSS DR7 ( Reid et "ar]|2010 ) in place 
of the BAO constraints, the full WMAP7 CMB power spectrum ( Larson et al.|[20TT ) in place of the 
shift parameters, and a prior on Hq from the SHOES (Supernovae and Hq for the Equation of State) 
program (Riess et al. 2009 2011). This Gaussian Hq prior, Hq = 73.8 ± 2.4 kms'i Mpc~\ makes 



use in its construction of many of the low-redshift z < 0.1 SNe la used in this paper. Their absolute 



magnitudes are calibrated directly using Cepheid variables in eight local SN la host galaxies ( Riess 



et al. 2011), the Cepheids themselves calibrated using different techniques: the geometric maser 



distance to the galaxy NGC 4258, trigonometric parallax distances for Milky Way Cepheids, and 
eclipsing binary distances for Cepheids in the Large Magellanic Cloud. The Hq prior was derived 



with SN la parameters from the MLCS2k2 distance estimator (Jha et al. 2007) given in Hicken 



et al. (2009bj). In principle, for complete consistency with this work, we would use the combined 
SALT2/SiFT0 fits for the same SN events to re-estimate Hq in a consistent way. However, the 



maser distance used in Riess et al. (2011) is not published, so we defer this exercise to a future 
analysis. However, we note that the uncertainty in Hq given by Riess et al. does include an 
allowance for the systematic of using SALT2 in place of MLCS2k2 (they quote an increase in Hq 
of 1.0 kms^^Mpc^^ with SALT2), so most of the systematic difference is likely already included 
in our error budget. 
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4. Results 

We begin by assessing the magnitude of the various systematic uncertainties in our analysis. 
For this, we use a simple cosmological model - a flat cosmology with a constant w - and the grid 
marginalization approach (§[3]). 

We then present our main cosmological results. We investigate a non-flat, w = —1 cosmology 
(fitting for Qm and fix), a flat, constant w cosmology (fltting for and w), a non-flat cosmology 
with w free (fitting for w, Qm and iljt), and a cosmology where w{a) is allowed to vary via a simple 



linear parameterization w{a) = wq + Wa{l — a) = wq + Waz/{1 + z) (e.g., Chevallier &: Polarski 



2001, Linder 2003), fitting for 0.^, wq, and Wa- We always fit for a, (3, and Mb- 



The confidence contours for and w ina flat universe can be found in Fig. [2] (upper left panel) 
for flts considering all systematic and statistical uncertainties. Fig. [2] also shows the statistical- 
uncertainty-only cosmological flts in the upper right panel. The best-fltting cosmological parameters 
and the nuisance parameters a, /3, and M%, for convenience converted to Mb assuming 

Hq = 70kms^^Mpc^^ (in the grid marginalization approach, Hq is not flt for as it is perfectly 
degenerate with Mb), are in Table [l] (for non-flat, w = —1 flts) and Table [2] (for flat, constant 
w flts). We also list the parameters obtained with the minimization approach for comparison. 
All the flts, with and without the inclusion of systematic errors, are consistent with a w = —1 
universe: we flnd w = — 1.043^qq55 (stat), and w = — LOGSlgggj (stat-|-sys). For comparison, 
with no external constraints (i.e., SNLS3 only) the equivalent values are w = — 0.90^q2o (stat) and 
w = -OmtoM (stat-Fsys) \Cll\. 



The lower right panel of Fig. [2] shows the importance of allowing the nuisance parameters a 
and (3 to vary in the flts, rather than holding them flxed at their best-fit values. This leads to 
not only smaller contours and hence underestimated parameter uncertainties, but also a significant 
bias in the best-fit parameters (Table [3|. Holding a and /? fixed gives w = — 1.117lQ'Qg2, ^ ~ O.60" 
shift in the value of w compared to the correct fit. 

The residuals from the best-fitting cosmology as a function of stretch and color can be found 
in Fig. [3j No significant remaining trends between stretch and Hubble residual are apparent, but 
there is some evidence for a small trend between SN la color and luminosity at C < 0.15 (indicating 
that these SNe prefer a smaller /3, or a shallower slope, than the global value). We examine this, 
and related issues, in more detail in § [5] 

Covariances between the nuisance parameters are small, with \r\ < 0.15 for most combinations 
of a, /3 and A4b- The exception is between Ai^^ and A^^, where the correlation is (as expected) 
larger (r ~ 0.6). Note that this positive covariance enhances the significance of the difference 
between and M'b beyond the simple statistical uncertainties listed in the tables. 
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4.1. Systematic error budget 



The ^m^w flat universe fits (Table[2]) represent a 5.2% statistical measurement of w and a 7.6% 
measure with systematics (i.e., ~ 5.5% with systematics only). The total systematic uncertainty 
is therefore comparable to, but slightly larger than, the statistical uncertainty. The full systematic 
uncertainty error budget can be found in Table [sj Systematic uncertainties generate about a ~70% 
increase in the size of the area of the SNLS3+BAO+WMAP7 ^m-w 68.3% confidence contour 
relative to a fit considering statistical errors only (compared to an 85% increase in the SN-only 



contour; see Cll ) 



The dominant systematic uncertainty is calibration (as in Cll| >, and in particular how well- 
known the colors and SED of the flux standard (BD 17° 4708) are - each of these two terms 
provides about a 20% increase in the contour area size over the statistical-only fit. The SNLS 
instrumental zeropoints and filter responses are also a large effect, generating a ~15% increase in 



each case. In part, this is because the SNLS data are calibrated to the Landolt (1992) system (for 



comparison to the low-redshift literature SNe), for which the color terms from the SNLS filters are 



large (Regnault et al. 2009). This situation should improve in the near future as new low redshift 



SN la samples observed in a similar filter system to SNLS become available, dramatically reducing 
these calibration uncertainties. 

By contrast, systematics caused by potential evolution in SN la properties (the parameters 



a and f5) are considerably smaller. As discussed in Cll, we find no evidence that a varies with 
redshift, and only marginal evidence for redshift variation in /3: explicitly da/dz = 0.021 it 0.07 
and d(3/dz = 0.588 it 0.40 (Cll). Is is unlikely that this (3 evolution is real (GIO); however, we 
conservatively adopt da/dz = 0.07 (the uncertainty in the slope) and d/3/dz = 1.0 in our systematics 
analysis. Even this amount of redshift evolution in a and /3 contributes an almost negligible effect 
(Table [s]) . The largest identified systematic uncertainty related to the astrophysics of SNe la is the 
implementation of the host-galaxy dependent term in eqn. 

The lower left panel of Fig. [2] shows the Vlm~w contours with all systematics included, except 
those related to calibration. These "no-calibration-systematics" contours are very similar to the 
statistical-only contours (only a factor 1.07 larger), with w = — 1.048ib]tQ Q5g. This represents a 
total error in w of ~ 5.5%, and a systematic contribution of ~ 1.8%, significantly smaller than when 
the calibration systematics are included. With our current knowledge and fitting techniques for 
SNe la, this represents the systematic floor given a negligible photometric calibration uncertainty. 



4.2. Cosmological results 

We now present our main cosmological results. We consider various combinations of the 
SNLS3, WMAP7, SDSS DR7 LRGs, and Hq datasets: WMAP7-hSNLS3+DR7 is the most similar 
to the constraints used in the grid marginalization approach (Table [2]) , but still differs as it uses 
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the full matter power spectrum of LRGs rather than the BAO constraint, and the full WMAP7 
power spectrum rather than the shift parameters. The best-fitting value of w (Table |4]) is therefore 
slightly different due to these differing external constraints even though the SN la constraints are 
identical, but the percentage error in w is the same at 7.6%. 



4-2.1. Constant w fits 

All the results are consistent with a spatially flat, w = —1 universe. Our results for a flat 
universe with a constant dark energy equation of state are 

ilm = 0.269 ±0.015 

_ _i Qci +0.069 

and, relaxing the assumption of spatial flatness, 

= 0.271 ±0.015 
= -0.002 ± 0.006 

W — i.UDy_QQg2, 

including external constraints from WMAP7 and SDSS DR7 and a prior on Hq (all quoted uncer- 
tainties in this section include both the SN statistical and systematic components). The confidence 
contours are in Fig. [4] and Fig. [5} and the corresponding best-fit cosmological parameters for various 
combinations of external constraints can be found in Table IH 

In Table [5] we give a full list of all the best-fit parameters from the constant w CosmoMC fits with 
the WMAP7, SDSS DR7 and Hq external dataset!|5 This includes some parameters which the SN 
la data do not directly constrain. and are the fractional energy densities of baryons and dark 
matter, r is the reionization optical depth, Ug is the scalar spectral index, Aq^ the amplitude of 
curvature perturbations at A; = 0.05Mpc~^, and ug the normalization of the matter power spectrum 
at 8h-^ Mpc. 

Of particular note is the high importance of the SN la dataset in placing meaningful constraints 
on w. Assuming a fiat universe, WMAP7±DR7 alone only measure w to ~20%, and adding the 
Hq prior (i.e., WMAP7+DR7±-ffo) only decreases this uncertainty to ~11%. Including the SNLS3 
dataset with WMAP7±DR7, by contrast, reduces the uncertainty to 7.7%. WMAP7+SNLS3 
together also provide a 7.7% measurement. With all external constraints, including the Hq prior, the 
measurement of w is 6.5%, comparable to WMAP7±SNLS3 alone. Note that the DR7 constraint 
has almost no effect on the uncertainty in the measurement of w - WMAP7±-ffo±SNLS3 has 
essentially the same uncertainty as WMAP7±DR7±iJo±SNLS3. 



^FuU parameter summaries for all combinations of external datasets can be found at https : //tspace . library . 
utoronto . ca/snls 
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The situation is slightly different when making no assumption about spatial flatness, but the ba- 
sic result of the high importance of the SN la data remains. In this case, the WMAP7+DR7+SNLS3 
data provide an 8.8% measurement, compared to 11.3% without DR7 (the DR7 data make impor- 
tant contributions towards constraining 0^). WMAP7+DR7+i?o alone can only make a ~ 20% 
measurement; adding SNLS3 improves this dramatically to ~ 8.5%. 



Variable w fits 

The final set of fits allow the equation of state parameter w to vary simply as a function of 
the scale factor, a, as tt;(a) = wq + Wail — a) , with a cosmological constant equivalent to = — 1; 
Wa = 0. We use a hard prior of wo+Wa < 0, from the constraint of matter domination in the early 
universe. The confidence contours are in Fig. [6] assuming a flat universe. The best-fit parameters 
are listed in Table [6j Again, we find no evidence of deviations from the cosmological constant. 
Assuming a flat universe we find 

^'m — 0.015 

WO = -0.905l°:J^6 



Wn 



-0.9841} 



1.094 
097- 



The Hq prior has only a small effect on the 'w{a) fits (Fig. [7]). As a comparison, with no Hq prior, 
we find Wo = — 0.9491q201) Wa = — 0.535l}']^'j'i- The SN la data are critical for a constraining 
measurement - without the SN data, fits for wq and Wa did not converge. 



4.3. Comparison to other results 



We compare our results to previous constraints on dark energy using SN la data. Komatsu 



the 



Hicken et al. 



et al. (|2011|), with a combination of WMAP7, BAO, the Riess et al. (|2009[) Hq measurement, plus 

Although 



(2009b) SN la dataset, found wq = -0.93 ± 0.13 and Wa 



-n 41+0-72 



this may appear to be better than our constraints, it did not include a proper handling of SN la 
systematics due to the lack of a consistent, published dataset at that time, and the uncertainties 
will therefore be under-estimated. Using a slightly larger SN la dataset ("Union2") the same 



WMAP7/BA0 constraints, and the |Riess et al.| ( |2009| ) Hq measurement, [Amanullah et aL] ( |2010 ) 
find = 0.2741°;°}^, O^, = -0.002 ± 0.007, and w = -1.052^°;°^^, comparable to our results 



(§ 4.2). However, those authors left a and /3 fixed when computing systematic uncertainties, which 
may underestimate the size of the final uncertainties (Cll, see also Fig. [2] and Table [3]), and likely 



under-estimated the magnitude of the photometric calibration uncertainties (see Cll for discussion). 



We also compare our results with those predicted by the Dark Energy Task Force (DETF; 



Albrecht et al. 2006) for an experiment of this type by calculating the figure of merit (FoM) for 
our combination of datasets. The exact definition of the FoM has some ambiguity: in the DETF 
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report it is defined as proportional to the reciprocal of the area of the error ellipse in the WQ-Wa 
plane that encloses 95% of the total probability, but the constant of proportionality is never stated. 



The values given in Albrecht et al. (2006) are based on performing a transform of variables from 



wq to Wp, the value of w at the so-called pivot redshift Zp, where Wp and Wa are uncorrelated. The 
FoM is then simply taken to be ^/ {(ywp'^wa)- This prescription is the most commonly used in the 
literature (e.g., Eisenstein et al.||2011 ), although some authors have been more literal in taking the 
area of the ellipse. FoM calculations should not assume a flat universe; however current data are 
poorly constraining without this constraint, so we follow the practice in the literature and assume 
flatness in our FoM numbers. 



For our WMAP7+DR7+SNLS3+ifo fit in a flat universe, we find Zr, ~ 0.19 and 



-1.063± 



0.082. Combined with our measurement of Wa, we find a FoM of 11.1 (see also Table [6]). Excluding 
the SHOES Hq prior gives a FoM of 10.6. Directly taking the reciprocal of the area of the 95% 
confidence intervals give 0.56 and 0.44 respectively. 

In DETF terminology, the SNLS3 sample represents a stage II SN experiment, and the final 
combination of this and other stage II experiments is predicted to give a FoM of ~ 50. However, 
these figures are difficult to compare with our results. The DETF calculations assume a far larger 
SN la sample of 1200 events (including, significantly, 500 at low-redshift), together with a (superior) 



CMB prior from the Planck satellite (Planck Collaboration et al. 2011) rather than WMAP, and 



also include cluster and weak lensing experiments. They do not include BAO or Hq information, 
although the latter only has a small effect on the FoM. 



5. Supernova sub-samples 

A complementary approach to checking and analyzing systematics, in particular those that 
are not susceptible to an analytical approach, is to break the SN la sample into sub-samples which 
probe either different experimental systematics, or different regions of parameter space of the SN la 
population. For example, the stretch-luminosity and color-luminosity relations are assumed to be 
linear, universal and invariant with (e.g.) SN properties, and the size of the SNLS3 sample allows 
us to test these assumptions in detail. Furthermore, the cosmological results should be robust to 
any segregation of the data if the systematics are handled correctly. 

For simplicity, we use the minimization approach in this section. We first split the SNLS 
SNe according to location on the sky (i.e., one of the four CFHT-LS deep fields). We then test 
the robustness of the nuisance parameters by splitting the sample by SN and their host galaxy 
parameters. 
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5.1. Segregation by SNLS field 

We first test for any variations in as a function of the field in which the SNLS SN occurred. 



SNLS observes in four fields distributed in right ascension (see Sullivan et al. 2006a, for the field 
coordinates), so we can test for a combination of photometric calibration systematics, as well as a 
physically varying w in different directions, by comparing the cosmological parameters we derive 
from each field. To do this, we adjust our minimization approach to fit for global nuisance 
parameters (a, /3, and A4 b) for all SNe, but with a different w and Vim in each of the four different 
SNLS fields. The average w and for the four fields is applied to the external SNe and for 
combination with WMAP7 and the BAO constraints. 

This approach effectively adds six new terms to the fit. The results can be found in Table [7| 
The only drops from 418.1 (for 466 DOF) to 414.9 (for 460 DOF), indicating that the fields 
are consistent. Comparing the individual field and w values to the average values of the four 
fields, and allowing for the covariances between the individual values, gives a of 3.38 for 6 DOF, 
consistent with the distribution and indicating no significant variability among the different 
fields. 



5.2. Segregation by SN properties 



The SN la light curve shape is well-known to vary systematically as a function of the SN 
environment. Fainter SNe la with faster light curves are preferentially located in older stellar 
populations, while the brighter examples with broad light curves tend to explode in late-type spiral 



or star- forming systems (e.g., Hamuy et al.|[l995 2000 Sullivan et al. 2006b). When coupled with 



the evidence for both a young and old component to the SN la progenitor population (Mannucci 



et al. 2005, 2006; Sullivan et al. 2006b Brandt et al. 2010), or at least a wide-range in the SN 



delay-times ( Pritchet et al.|[2008 Totani et al.||2008 ), a natural prediction is a subtle change in the 



mix of SN light curve shapes with redshift ( Howell et al.|[2007 ). 



If the stretch-luminosity relation is universal across SN stretch and progenitor age, this pre- 
dicted drift will not impact the determination of the cosmological parameters - a low-stretch SN 
should correct equally well as a high-stretch SN. To test this, we split our SN sample into two 
groups (s<l and s>l) and perform independent cosmological fits to each sub-sample. This obvi- 
ously restricts the lever-arm in stretch and so the a coefficient is less well-determined; nonetheless 
this is a useful test of the utility of SNe la across different environments. 

When considering sub-samples of SNe classified by SN properties (for example stretch or color) , 
the Malmquist corrections that are applied globally to the sample will not be appropriate. For 
example, low stretch SNe la are intrinsically fainter and will suffer from a larger selection effect 
at high-redshift than high stretch events. Rather than apply different Malmquist corrections for 
these different sub-samples (which could in principle be derived from simulations such as those 
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in Perrett et al. 2010), we instead restrict the SNe to the redshift ranges over which selection 
effects are reduced. For the SNLS sample, we restrict to z < 0.75, and for the SDSS sample, we 
restrict to z < 0.3. These data also tend to be better observed, being brighter, with smaller error 
bars on the SN parameters. We also discard the small HST sample. We perform both statistical 
and statistical+systematic uncertainty fits - the comparison is useful as we are interested in the 
differences between SN la sub-samples, and many of systematics affect sub-samples in a similar 
way. 

The results are given in Table [8] and Fig. [8j the latter generated from the covariance matrices 
of the fits.. The derived cosmological parameters are consistent between low and high-s SNe la, 
although the nuisance parameters show some differences: the a values are consistent between low-s 
and high-s, but much larger than the full sample. The Mb values show large differences between 
low-s and high-s. However, there are significant (and expected) covariances between a and Mb in 
these fits: for the low-s group, Mg increases as a decreases, while for the high-s group M^ decreases 
for decreasing a (and a similar trend is seen for M^). As Mb is defined at s = 1, which the SNe 
in neither group sample, the two parameters become quite inter-dependent. When considering the 
joint confidence contours between a, M^, and (Fig. [8| only mild tensions are seen. A slightly 
larger tension is seen in the value of /? for the two sub-groups: the values 2.90 it 0.16 and 3.25 it 0.14 
differ at ~ 2.2c7 (Table |8j). 

In a similar vein, we also split the sample by SN color at C = 0. Unlike stretch, SN color 
(examined independently of SN luminosity) shows no clear evidence for variation with environment. 
No robust trends have been found despite rigorous examinations (e.g., Sullivan et al.|2010 Lampeitl 



et al.pOlOD , and we also find no significant trends of the nuisance parameters varying as a function 



of SN color. 



5.3. Segregation by host galaxy characteristics 

We also segregate the SNLS3 sample according to the environment in which the SN exploded. 
This is a fundamentally different test, dividing the sample according to the environment of the SN 
rather than its direct properties (although significant correlations exist between, for example, SN 
stretch and host galaxy Mgteiiar or star- formation rate). For these tests we use Mgteiiar to segregate 



the sample, as used in the cosmological fits, following Sullivan et al. (2010). As we already use 
a different Mb for SNe la in low and high Mgteiiar galaxies in our standard fits, we simply adapt 
this approach to additionally fit for a different a and /3 in the two host classes. This has the 
advantage that we are no longer comparing the results of different fits as we fit for all the nuisance 
parameters simultaneously, and we can use the full systematics covariance matrix, as well as the 
full SN la sample. (We also provide the fit results when we physically divide the sample into two, 
for comparison with § |5.2l in Table pi) Note that in these fits we do not include the host galaxy 



systematic term, as we are trying to examine the effect of any host galaxy dependence. 
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The "multi nuisance parameter" fit results are shown in Table [9j which gives the values of the 
nuisance parameters themselves, and Table 10, which gives the effect on w and the x^. The joint 



confidence contours for some of the nuisance parameters are shown in Fig. [9j For completeness, we 
also give results when only one Mb is used. 

The data do not support the addition of a different a parameter in low and high mass galaxies 
- when this is added to the fits, the a values are generally consistent and the quality of the fit, 
as indicated by the x^, is unchanged. This is true even when only one Mb is used in the fits, 
and suggests that a is fairly insensitive to the details of the environment and characteristics of the 
SN la progenitor stellar population. 

However, there is evidence for different Mb (as already fit for) and different (3. The value of f3 
is ~ 3.7 in low mass galaxies, versus ~ 2.8 in high-mass galaxies, regardless of whether two A4b are 
used. A similar trend is seen when physically dividing the sample into two and performing separate 
independent fits as in the previous section (Table [8|, and is consistent with the /3 difference seen 
between low and high-stretch SNe in § |5.2[ as low-stretch SNe are preferentially found in massive 
host galaxies. Generally, the two /3 values show only a very small positive covariance, and differ 
at the ~ 4.3cr level. There is also a substantial reduction in the of the fit when including two 
/3 terms. For example, fitting for two /3s and two Mb reduces the to 405.4 from the 423.1 
obtained if only one f3 is used (for 465 and 466 DOF). An F-test indicates that this additional term 
is required at ~ 4.5(T. 



The variation of f3 with host properties has been observed at lower significance by [Sullivan 



et al. (2010) in the SNLS sample, and by Lampeitl et al. (2010) in the SDSS SN la sample (a larger 



SDSS SN la sample than the one used in this paper). Lampeitl et al. (2010) find (3 ~ 2.5 in passive 
host galaxies, and (3 ~ 3.1 in star-forming host galaxies. Using their full sample, the significance 
is ~ 3.5fT, although this drops to < 2a when considering only cosmologically useful events with 
normal stretches and colors using similar cuts to those used in this paper. 



5.4. Discussion 

The most significant result from our analysis of the SN sub-samples is the additional variation 
of (3, as well as Mb, between low- and high-Msteiiar host galaxies. This effect appears real in 
our data and so should be accounted for appropriately in our cosmological results. We therefore 
examine the systematic effect of not including this term, and compare to our existing systematic 
uncertainty error budget. 

Compared to a statistical uncertainty only fit, the addition of two as and two /3s, gives a 



A^lm = 0.00 and Aw = 0.005 (Table 10). The mean statistical-only errors on and w are 0.0148 
and 0.0545 (Table [s]). Adding in quadrature the shifts measured when including the two as and /3s 
increases the w uncertainty to 0.0547. This total uncertainty on w is smaller than the uncertainty 
obtained when including host galaxy systematic term listed in (0.0559; Table [s]), i.e. the effect on 
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w of introducing two /3s is smaller than our current host galaxy systematic term. 



Note this would not be the case if we had neglected all nuisance parameter variation, i.e. had 
only used one Mb in our cosmological fits. In this case, AQm = 0.005 and Aw = 0.055 (the 
shift from considering one Mb to considering two Mb)', in the case of w this is a shift larger than 
our statistical uncertainty, and is comparable to our total systematic uncertainty, becoming the 
dominant term in the error budget. Thus while the use of different Mb is essential for a SN la 
cosmological analysis, the use of two /3s and two Mbs is not. Note that no previous SN la cosmo- 
logical analysis has performed this host galaxy correction, indicating that systematic uncertainties 
will be significantly under-estimated in these studies. 

A similar argument can be made using the fits including systematic uncertainties. Here, we 



compare fits that do include the host galaxy systematic term (unlike the numbers in Table 10), 
as we wish to examine whether the size of any shift in the cosmological parameters with the 
addition of two /3s is accounted for by our existing systematic uncertainty error budget. In this 
case. Aw = 0.015. Our total error in w is 0.0810, compared to 0.0800 excluding the host systematic 
term (see Table [s]). Adding the 0.015 in quadrature to this 0.0800 gives a w uncertainty of 0.0813, 
a total uncertainty on w almost the same as that obtained when using the host systematic term. 
Therefore, our conclusion is that while the two /3 effect appears real in our data, it is adequately 
accounted for by our systematic uncertainty error budget. 

Although the variation of (3 with host parameters is not a concern for this cosmological analysis, 
it does have implications for the physical origin of color variation in SNe la which may impact future 
surveys. A long-standing observation is that the slope of the relation between Mb and B — V (i.e., 
13) is <C 4.1 ( |Tripp||1998 | [Astier et al^|2006D , the value expected based on Milky Way like dust if /3 
is interpreted as the ratio of total-to-selective extinction Rb (where Rb = Rv + 1; and Ry ~ 3.1 
for the Milky Way). The effective j3 for SNe la is likely a confiation of different physical effects. 



including extinction by dust (which may vary with host type, e.g., Lampeitl et al. 2010 Sullivan 



et al. 2010, and this paper) and intrinsic variation in SN color that does not correlate with SN 
light curve shape (e.g., Folatelli et al.|[2010 ), and which may depend on variables such as explosion 



asymmetry or observational viewing angle (Kasen et al. 2009 Maeda et al. 2011b). Recent work 



has shown that SN color is also correlated with SN la spectral features, with SNe possessing faster 
ejecta velocities having redder colors at fixed Mb (or equivalently brighter Mb at fixed color) in 
samples with very red SNe excluded ( Foley fc KasenpOll ). 

Under the assumption that any intrinsic SN la color-luminosity relation has a smaller effective 
j3 than that from dust (as seems likely given we observe (3 < 4.1), our results are qualitatively 
consistent with a scenario in which dust extinction modifies this intrinsic color-luminosity relation. 
The lowest Mgteiiar host galaxies are those with the largest specific SFRs, and therefore the largest 
dust content. We would therefore expect to find SNe la with a larger effective /3 (i.e., closer to the 
true dust value) in lower Mgteiiar hosts, which is consistent with our observations. In more massive, 
passive host galaxies, we are likely observing a /3 closer to the intrinsic value. 
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6. Conclusions 

In this paper we have presented the cosmological results for the Supernova Legacy Survey 



(SNLS) three-year SN la sample (SNLS3; GIO Cll), combined with other constraints from the 
literature. Our SN la sample contains 472 SNe, including 242 from SNLS (GIO), 123 at low 



redshift, 93 from SDSS-SN, and 14 from HST (see Cll). We have performed analyses investigating 
the cosmological parameters Qm, ^fcj and the dark energy equation-of-state parameter w. A key 
aspect of our analysis is the inclusion of all identified SN la systematic uncertainties directly in 



our cosmological fits (Cll). The inclusion of these systematic uncertainties has two key effects. 



Foremost, the uncertainties that we quote on the cosmological parameters reflect the systematic 
component. Second, correlations in brightness, stretch and color between different SNe due to their 
being affected by the same systematic are accounted for during the cosmological fitting stage. We 
also correct for recently identified trends between SN la brightness and host galaxy stellar mass, 
and account for the effect of systematic differences from the use of two independent SN la light 
curve fitters. 



Our main results are: 



1. For simple cosmological fits assuming a flat Universe and constant w, combining the SNLS3 
sample with BAO observations and the WMAP7 CMB "shift" parameters gives = 
0.276]']Q'gJ3 and w = —I.OAS'^q'q^'^, where the error is statistical only. When we include 
ah identified SN la systematics in the fits, we find Qm = 0.274+°;°J^ and w = -1.068+|5;[J^^. 

2. In terms of the contribution towards the uncertainty in measuring w in the above fits, our 
systematic and statistical uncertainties are approximately equal (5.5% and 5.2% respectively). 
However, the systematic uncertainty error budget is dominated by the photometric calibration 
of the SN fluxes, rather than uncertainties related to the astrophysics of the SNe themselves. 
Neglecting calibration uncertainties, likely to be dramatically reduced in the future, gives a 
systematic uncertainty of ~ 2%. 

3. When including the SHOES prior on Hq, together with the full WMAP7 power spectrum and 
the power spectrum of LRGs in SDSS DR7, we find = 0.269± 0.015 and w = -IMl^S 
using the CosmoMC fitter, a 6.5% measure of w. When we relax the assumption of a flat 
Universe, we find Qm = 0.271 ± 0.015, Qk = -0.002 ± 0.006, and w = -1.069+[]:|]^2- These 
include all SN systematic uncertainties. 

4. We consider a simple parameterization of the time variation of w as w{a) = wq + Wa{l — a). 
Assuming a flat Universe we find = 0.27ll[!:[}i^, wq = -0.905l[]:}^J^ and Wa = -0.9841};°^^. 
This includes WMAP7, SDSS-DR7, and the SHOES Hq prior as external constraints. Our 
results are equivalent to a DETF figure-of- merit of ~ 11. 



5. We investigate astrophysical systematics in our SN la sample by breaking it into sub-samples 
based on SN light curve and host galaxy parameters. Cosmologies determined from these SN 
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sub-samples are fully consistent. However, we find significant evidence (4.4(t) for a different /3 
between low and high stellar mass host galaxies (as well as a different Mb, which we already 
account for) . The effect of this varying (3 on the cosmology lies well within the current system- 
atic uncertainty assigned to host-galaxy-dependent corrections, but our analysis emphasizes 
the critical need to make host galaxy related corrections when determining the cosmological 
parameters. 

When the SNLS3 sample is combined with measurements of large scale structure, observations 
of the CMB, and a prior on Hq, the constraints on dark energy presented here are the tightest 
available, and directly include all identified SN la systematic uncertainties in the analysis. All our 
results are consistent with a flat, w = —1 universe. 

The primary contributor to the systematic error budget is the calibration of the SN la fluxes, 
both placing them on a consistent system between different SN la surveys, and then interpret- 
ing that system when fitting the SN light curves. The magnitude of the identified astrophysical 
systematics is significantly smaller than calibration related uncertainties. Those that have been 
uncovered in the SN la population can be adequately controlled using empirical corrections based 
on the properties of the SN la host galaxies. 

The implication is that, if the calibration-related systematics can be reduced, SNe la are a 
long way from being systematics limited. In part, the calibration-related systematics arise from the 
need to calibrate the griz SNLS and SDSS filter sets to the UBVR system used for the majority of 
the low-redshift SN la data. This situation is set to improve considerably as improved low redshift 

. These new samples will be 
directly calibrated to the SNLS system (or vice-versa) and will eliminate, or at least substantially 
reduce, the main systematic uncertainties, allowing the full potential of the SNLS sample to be 
unlocked. 
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Fig. 1. — SNe la from the SNLS3 sample in color-color space. The colors are measured at maximum- 
light in the rest-frame i?-band. SNe are coded according to the sub-sample from which they are 
drawn with the mean redshifts shown, and the best-fitting SiFTO color-color law is over-plotted 
for each sample - see Conley et al. (2008) for details of this fit. The arrow indicates the direction 



of the Milky Way extinction vector using a Cardelli et al. (1989) law. The CfA sample shows an 
offset from, and is statistically inconsistent with, the other samples. 
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Fig. 2. — Confidence contours in the cosmological parameters and w arising from fits to the 
combined SN la sample using the marginahzation fitting approach, illustrating various systematic 
effects in the cosmological fits. In all panels, the SNLS3 SN la contours are shown in blue, and 
combined BA0/WMAP7 constraints (Percival et al. 2010 Komatsu et al. 2011) in green. The 



combined constraints are shown in grey. The contours enclose 68.3%, 95.4% and 99.7% of the 
probability, and the horizontal line shows the value of the cosmological constant, w = —1. Upper 
left: The baseline fit, where the SNLS3 contours include statistical and all identified systematic 
uncertainties. Upper right: The filled SNLS3 contours include statistical uncertainties only; the 
dotted open contours refer to the baseline fit with all systematics included. Lower left: The filled 
SNLS3 contours exclude the SN la systematic uncertainties related to calibration. Lower right: 
The filled SNLS3 contours result from fixing a and /? in the cosmological fits. See Table [2] and 
Table [3] for numerical data. 



0.7 



0.9 1.0 
SN Stretch (s) 



-0.2 



-0.1 



0.0 

SN Color ic) 



0.1 



0.2 



Fig. 3. — Residuals (in magnitudes) from the best-fitting flat cosmology as a function of stretch 
(left) and color (right). Residuals are defined as niB — m'^°'^, i.e., negative residuals indicate 
brighter SNe (after application of stretch- and color-luminosity relations). Red points show the 
mean residuals in bins of stretch and color. The dashed line indicates a zero residual. 



- 33 - 



WMAP7 + ... 




0.1 0.2 0.3 0.4 0.5 



Fig. 4. — Confidence contours in the cosmological parameters 0,m and w assuming a flat universe, 
produced using the CosmoMC program. The SNLS3 contours are in blue, the SDSS DR7 LRG 
contours in green, and the Hq prior in red. WMAP7 constraints are included in all contours. The 
contours enclose 68.3% and 95.4% of the probability and include all SN systematic uncertainties. 
The dashed line indicates w = —1. Numerical results are in Table [IJ 
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Fig. 5. — Confidence contours in the cosmological parameters f^m, ^de^ ^fci and w produced using 
the CosmoMC program. The SNLS3 contours are in blue, the SDSS DR7 LRG contours in green, 
and the Ho prior in red. WMAP7 constraints are included in all contours. The contours enclose 
68.3% and 95.4% of the probability and include all SN systematic uncertainties. In the left-hand 
panel, the dashed line indicated w = —1; in the center and right-hand panels the line indicates a 
flat (r^fc = 0) universe. Numerical results are in Table ll 



SNLS3 + SDSS DR7 LRGs+WMAP7+i7n (Flat) 




Fig. 6. — Combined confidence contours in Qm, wq, and Wa using SNLS3, WMAP7, SDSS DR7 
LRGs, and a prior on Hq. A flat universe is assumed, and we enforce a prior oi wq + Wa < - any 
apparent discrepancy with this prior is a result of smoothing the CosmoMC output. The horizontal 
dashed lines indicate a cosmological constant (wq = —1; left) and a non-varying w {wa = 0; right). 
All SN la systematic uncertainties are included. Numerical results are in Table [6j 
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Fig. 7. — The effect of the Hq prior on the wq, and Wa contours in a flat universe. The red 
contours show the fits with the Hq prior, and the blue contours without. 
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Fig. 8. — Joint confidence contours between the nuisance parameters a, (left), and (right) 
for low-stretch (solid line) and high-stretch (dashed line) SNe la, using the restricted SN la sample 



described in § 5.2, The contours enclose 68.3% and 95.4% of the probability, and the fits include 
aU systematic uncertainties. Only mild tensions exist. 
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Fig. 9. — Joint confidence contours in the nuisance parameters /3 and Mb (left) and a and 
Mb (right), allowing all the nuisance parameters to vary according to host galaxy stellar mass. 
ai//3i/M^ refer to SNe la in hosts with Msteiiar < lO^^M©, and a^j^ilM^ to SNe la in hosts with 
Msteiiar > 1O^'^M0. The full SNLS3 sample is used, and all systematic uncertainties are included. 
A significant variation in /3 with host Msteiiar is observed. 
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Table 1. Cosmological results assuming w = —1 for the SNLS3 sample plus BAO^ and WMAP7^ 



Fit 


a 




M}^^ Ml"" 








Marginalization fits 


Stat only 
Stat + sys 


i.40i_o.lOO 
1 498+0-121 


Q lfi^+0-105 

a.lUd_Q Qgg 

q 9^0+0.121 


-19.122 -19.198 
-19.144 -19.196 


n 97r;+0.016 

n 97q+0-0l9 
^•^'^-0.015 


n 797+0.015 
'^■'^'-0.013 

Q 724+0-017 

'^^-0.016 


minimization fits 


Stat only 
Stat + sys 


1 Qaq+0.085 
1 Qf;c+0.086 

i-'jUO_Q Qg4 


o 144+0.095 
q 1 09+0.102 

a.lOZ^Q Qgg 


-19 121+0-015 _iq iqfi+0.013 
-19 162+° °28 _iq 2n6+°-°24 

iy.iDZ_QQ2g iy.ZUD_QQ24 


n 97q+0.015 

U.Z(0_QQ;^4 

n 274+0-017 

"-'•^'^-0.016 


0.7291°:°}^ 

n 7q9+0.016 



,-1 



'^For an s = 1 and C = SN la. Computed from 7W_b (§ [3j) assuming Hq = 70kms^ Mpc" 
Errors on M b are not available in the marginalization (grid) approach as the variable is analytically 
marginalized; the quoted value is an estimate only. 



'^Percival et al. 



(2010) 



^Using the WMAP7 "shift" parameter R, the "acoustic scale" la-, and the decoupling redshift 



from Komatsu et al. (2011) 
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Table 2. Cosmological results assuming a flat universe and constant w for the SNLS3 sample 

plus BAO and WMAP7 



Fit 


a"" 




Mh Ml 






W 






Marginalization fits 


Stat only 
Stat + sys 


1 Q«7+0.086 

i.OO/_QQg4 


o 1^4+0.096 

+0.101 
"^•^ ' ^-0.099 


-19.164 -19.227 
-19.175 -19.220 




276+°-°^^ 

+0.019 
'-'•^'^-0.015 


-1 043+°-°^^ 

1 nfi8+0-080 
-i.0b»_oog2 


minimization fits 


Stat only 
Stat + sys 


1 QQc;+0.085 
1 QR7+0.086 

i.OU/_QQg4 


q 1/IO+0.095 

+0.101 
"^•^ ' ^-0.099 


-19.130j:[!:j]i^ -19.203j 
-19.155l[!:[]^^ -19.200j 


-0.016 
0.016 

-0.023 
0.023 


"-'•^'^-0.014 

r, 979+0.017 
'-'•^'^-0.016 


-1.039 
-1.058 


f0.052 
-0.055 

f0.078 
-0.082 



'^Note that the values of these nuisance parameters differ very slightly from those from the SN-only 
fits given in Cll due to small correlations between the cosmological parameters and the nuisance 
parameters. 
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Table 3. Detailed summary of systematic uncertainties 



Source 



w 



Relative area*^ 



Statistical only 


n O7fiq+O.oi63 
^•^'""^-0.0132 


-1 n43n+°°5^^ 


1.0 


All systematics 


n 97qfi+0.0186 


i.UDm_Q Q821 


1.693 


All systematics, except calibration 


0.2756lHil^ 




1.068 


All systematics, except host term 


0.2738inil^ 


i.UD44_Q Q809 


1.677 


All systematics, fixing a, 


0.2656lHiII 


1 11^8+00807 


1.641 



Contribution of different systematics: 



Calibration 




_ _ _ ^ _ _i_n n'V'VA 


1.614 


SN la model 


0.2767l°:gii 


-1.04031°:°^^^ 


1.013 


Peculiar velocities 


1 r\ ni art 

2761+°-°^^^ 


— 1 04^9+00544 
i.U^OZ_o 0548 


1.002 


Malmquist bias 


n ni-rr'n+n 0163 

0.2758l{]:^i^^ 


-, r\ .„ . 4-0 0548 

-1.0474lJ]:^548 


1.014 


Non SN la contamination 


0.2763+Hi3i 


-1.0430+H^4^ 

— U.U54D 


1.000 


Milky Way extinction 


0.2762+0-0J64 


-1.0441+0-|]i3 


1.023 


SN redshift evolution 


fi 97f?q+0.0163 


-1.0408tn^r7 


1.017 


Host galaxy term 


n O7f:o+0-0163 
^■^'"^-0.0132 


-1 n4^3+°°556 

i.U^OO_o 0562 


1.029 


Calibration: 








Colors of BD 17° 4708 


271 Q+0 0170 


1 0790+0-0639 
l.U/ZU_o 0639 


1.239 


SED of BD 17° 4708 


r. 977-1+0.0170 
u.z/ / -L_o.oi38 


±.uo»u_o 0630 


1.205 


SNLS zeropoints 


n O7f;7+0.0168 


-1 0421+°°^°^ 


1.166 


Low-z zeropoints 


n O7c;q+0.0164 


1 nf^97+0-0578 
l.UOZ(_o 0586 


1.078 


SDSS zeropoints 


-0.0133 


-1 0411+°-°^'^^ 

l.U^ll_0 0548 


1.015 


SNLS filters 




1 n'?'?n+0-0585 

i.UOOU_o.0586 


1.136 


Lowz filters 


U.zmO_QQ^32 


i.U4UZ_o 0550 


1.010 


SDSS filters 


+0.0164 
U.Z( /U_Q 0133 


i.UOyD_o.05.18 


1.007 


HST zeropoints 


u.z (u»_o 0133 


-1 041 2+°-°5'^'* 

i.U+iZ_o 0548 


1.007 


NICMOS nonlinearity 


2767+°-°^^^ 

U.Z/D^_Q Q133 


1 0418+0-0545 
-i.U4i8_oo548 


1.009 


SN la model (light curve fitter): 








SALT2 vs. SiFTO 


n O7f;7+0-0163 
U.Z/D^_Q 0132 


-1 0404+°-°^'^3 

l.U^U4_0 0547 


1.012 


Color uncert. model 


fi 97^0+0.0163 
U.Z /oo_o 0132 


-1 0430+°°^^3 

l.U^OU_o 0546 


1.001 



SN la redshift evolution: 



a 



n 97«q+0.0163 
u.z /uo_Q 0132 



-1 0430+°°^^^ 
i.U^OU_o 0546 



1.000 
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Table 3 — Continued 



Source Qm 


w 


Relative area*^ 


(3 0.2763tlS 




1.017 



'^The area of the Q,m~w 68.3% confidence contour relative to a fit considering statistical errors 
only. The contours are computed with the marginalization (grid) approach, include BAO and 
WMAP7 constraints, and assume a flat universe; see text for details. 

^All our cosmological results have a and /3 free in the fits. This entry shows the effect of 
incorrectly holding them fixed. 



Table 4. Cosmological results'' obtained using the CosmoMC fitter with a constant dark energy 

equation of state 



Parameter WMAP7+SNLS3 WMAP7+DR7 



WMAPT+i/O*^ WMAP7+DR7 



WMAP7+DR7 WMAP7+Ho 
+SNLS3 +SNLS3 



WMAP7+DR7 
+Ho+SNLS3 



Flat, constant w. 

W 



0.262 



+0.023 



-1.016 



+0.077 



Ho 

Non-flat, constant w 



71.58l|;42 



W 

Ho 



+0.050 



0.259 
001+° ''i5 

-l.Ul»_o.iio 



0.329 
-0.826 
64.42 



+0.034 
0.039 
+0.166 



+4.25 
4.38 



0.312 
-0.006 
-1.027 

66.36 



+0.051 
0.051 

+0.013 
0.012 

0.379 
0.386 
+5.73 
-5.87 



246+° °20 

-1 114+0'i" 
J^.ll't_0.113 

74 1 1 +2-58 



267+° °l'^ 
'^■^0' -0.017 

-1 110+°-l22 
^■^^^-0.120 

79 01 +2.46 
/ z.zi_2 42 



„ 91-0+0.021 

U.ZOO_Q Q21 

-0 012+°"°* 

-1 445+0-298 
^■^^■^-0.292 

73.54l|j^ 



0.2843 



-1.021 
69.77 



+0.078 



+2.07 
2.07 



U.za4_Q Q21 
-0 009+°-°°** 

-1.0b8_Q pgg 

67.85l|j? 



250+°-°" 
_1 n'^7+0-068 

l.UO(_Q Qgg 

72.85j;i;!^? 



r, 247+0.018 
'J.^*' -0.018 

004+°-°°'^ 

-1 001+°-°^^ 
l.UUi_(j 032 



73.64 



+2.24 



269+° °^5 
-1 061+°°'^'' 

l.UDl_Q ggg 

71 57+1-65 
'-'-•3'-1.65 



r, gyi +0.015 

-0 002+° °°^ 
-1.069lg:°g^i 
7l.l8li;f7 



''The values quoted are the expectation values of the marginalized distributions, not the best fits, with the 68.3% marginalized 
values quoted as the errors. All SN systematic uncertainties are included. Note that as the non-SN constraints used in 
CosmoMC differ slightly from those used in Table [5] the cosmological parameters are different. The closest comparison is 
SNLS3-l-WMAP7-fDR7. 



''We show only flat universe fits for the WMAP-|-_ffo combination; the fits were not constraining for non-fiat cosmologies, 
with the lower bound on w unconstrained. 
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Table 5. The full set of cosmological parameters obtained with the CosmoMC fitter 



Class 


Parameter 


Const, w 


Const, w 






flat 


non-flat 



Primary 




r> r>ro +0 054 


„ r»/^r4-0 056 

^•^"^-0.056 




Qc 

a 
(7 


nil ^n+0.0041 

0-1149^^.0041 
i.Udy»_ooo27 


-, -, A r+0.0047 

0-1145Io.oo47 
1 nzLm +0-0026 




T 


0.0871^:^7 


0.0881°:^°? 




Qf, 




U.UUZ_Q QQg 




Wo 


-1 061+°°^'' 

l.UD±_Q Qgg 


l.UU»_Q Qg2 




Us 








log[10lOAo5] 


q flntc+0.033 


004+0.033 




a 


i-45il°:l^^ 






P 




ocrq+0.111 


Derived 




n 701+0.015 

'•5J--0.015 


r. 70-1 +0.015 

'-'•'"^-'^-0.015 




Age^ 


13.7ll°:}J Gyr 


13.781°:^! Gyr 






0.2691°:°}^ 


n 071+0.015 

-'--0.015 




0-8 


0.8501^1 


0.847l°:°i 




^re 


10.5511:?° 


10.5511:?° 




^0 


71.571^:^1 km s-i Mpc-^ 


71.18ll:p km s-i Mpc-i 



^The current age of the universe 

'^The redshift at which the reionization fraction is a half 
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Table 6. Cosmological results obtained with CosmoMC assuming a variable dark energy equation 

of state and a flat universe 



Parameter 



WMAP7+DR7 
+SNLS3 



WMAP7+DR7 
+SNLS3+i7o 
(stat. only) 



WMAP7+DR7 

+SNLS3+i/o 
(stat+sys) 



Wo 
Wa 

Ho 

FoM^ 



n 000+0.019 
'-'•^"^-0.019 

-n 949+0-198 
u.y^y_o.2oi 



-0.535 



-1.109 
-1.111 



70.26111° 
10.6 



274+°-0i^ 

'^■^'^-0.014 



-0.870 



-1-0.139 
-0.139 



-0.9381^21 

71.381}J° 
21.5 



0.271 



-fO.015 
-0.015 



-0.905l°:Ji 
-0.984l};°^^7 

11.1 



The DETF (Albrecht et al. 2006) figure of merit (FoM), implemented here as l/((T^pCr^„); see 



4.3 



Table 7. Cosmological results fitting for a different w/flrn in each SNLS field 



Fit 




w 


a 






r.m.s. 


Basic fit 


n 979+0.017 
U.z ( ^_o.oi6 




1 qfi7+0.086 
1.00(_o 


o iyq+0.102 
O.i ( »_o.o99 


418.1 


0.153 


4 field fit^ 


U.ZDO_Q^33, 
n oil +0.095 
U.O±l_Q ]^23' 

0.2681°:!^° 


l.U^'±_g 300, 

1 9qr+0.312 

n Q1^ +0.222 

u.yoi_g 245, 
l.UOO_g 310 


1 QC,R+0.086 
i.oou_o 084 


q 1 oq+0.103 
O.i00_o 101 


414.9 


0.153 



'^Fitting for global nuisance parameters and a different and w in each of the four SNLS fields, 
average values of the four il„i and w values are applied to non-SNLS SNe. See § |5.1|for details. 



The 
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Table 9. Nuisance parameter variation for low and high Mgteiiar host galaxies. The a columns 
give the significance of the difference in that nuisance parameter between the two host types. 



Fit 


q1 




CTdiff- 










0"diff 








Ml 




Statistical: 




























la,ip,lMB 


1.29 ±0.08 






3.15 ±0. 


,09 








-19. 


,188 ±0.015 








la,l/3,2Ms 


1.39 ±0.08 






3.14 ±0. 


,09 








-19. 


,130 ±0.019 


-19, 


,203 ±0.016 


5.0 


la,2/3,lMs 


1.27 ±0.08 






3.70 ±0. 


,16 


2.78 ±0. 


11 


4.7 


-19. 


,186 ±0.015 








2q,1/3,1Mb 


1.10 ±0.14 


1.34 ±0.09 


1.5 


3.16 ±0. 


,09 








-19. 


,191 ±0.015 








la,2/3,2Ms 


1.37 ±0.08 






3.64 ±0. 


,16 


2.81 ±0. 


11 


4.3 


-19. 


,130 ±0.019 


-19, 


,197 ±0.015 


4.7 


2a,l/3,2Ms 


1.43 ±0.17 


1.39 ±0.09 


0.2 


3.14 ±0. 


,09 








-19. 


,128 ± 0.021 


-19, 


.203 ±0.016 


4.5 


2a,2/3,lMB 


1.00 ±0.15 


1.34 ±0.09 


2.1 


3.74 ±0. 


17 


2.79 ±0. 


11 


5.0 


-19. 


,190 ±0.015 








2a,2/3,2Ms 


1.33 ±0.17 


1.38 ±0.09 


0.3 


3.64 ±0. 


16 


2.81 ±0. 


11 


4.3 


-19. 


,133 ± 0.021 


-19, 


,198 ±0.015 


4.0 


Statistical+Systematic^ : 


























Ia,l/3,1MB 


1.28 ±0.08 






3.19 ±0. 


10 








-19. 


,193 ± 0.022 








la,l/3,2MB 


1.39 ±0.09 






3.18 ±0. 


10 








-19. 


,130 ±0.025 


-19, 


,206 ±0.022 


4.6 


la,2^,lMs 


1.27 ±0.08 






3.75 ±0. 


,17 


2.80 ±0. 


12 


4.5 


-19. 


,179 ± 0.021 








2a,ll3,lMB 


1.08 ±0.14 


1.33 ±0.09 


1.6 


3.21 ± 0. 


10 








-19. 


,196 ± 0.022 








la,2/3,2Ms 


1.36 ±0.08 






3.65 ±0. 


17 


2.85 ±0. 


12 


4.0 


-19. 


,127 ± 0.024 


-19, 


,198 ±0.021 


4.6 


2a,l/3,2Ms 


1.41 ±0.17 


1.38 ±0.09 


0.2 


3.18 ±0. 


10 








-19. 


,129 ± 0.026 


-19, 


,206 ±0.022 


4.6 


2a,2/3,lMs 


0.98 ±0.16 


1.33 ±0.09 


2.1 


3.78 ±0. 


17 


2.81 ±0. 


12 


4.7 


-19. 


,185 ±0.021 








2a,2/3,2Ms 


1.31 ±0.18 


1.37 ±0.09 


0.3 


3.65 ±0. 


,17 


2.85 ±0. 


12 


4.1 


-19. 


,129 ± 0.025 


-19, 


,198 ±0.021 


4.2 



^We exclude the host galaxy systematic term in these fits; see § |5.3| for details. 
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Table 10. Cosmological fits including multiple nuisance parameters for low and high Msteiiar host 

galaxies 



Fit w r.m.s. xVDOF 



Statistical: 



la,ip,lMB 


0.269 ± 


014 


-1 


094 







055 





153 


453 


6/467 


la,l/3,2MB 


0.274 ± 


014 


-1 


039 







054 





152 


429 


6/466 


la,2P,lMB 


0.271 ± 


014 


-1 


072 







054 





152 


430 


2/466 


2a,l/3,lMB 


0.270 ± 


014 


-1 


088 







055 





152 


450 


8/466 


la,2P,2MB 


0.274 ± 


014 


-1 


034 







053 





150 


410 


9/465 


2a,lf3,2MB 


0.274 ± 


014 


-1 


039 







054 





153 


429 


7/465 


2a,2f3,lMB 


0.272 ± 


014 


-1 


064 







054 





151 


425 


7/465 


2a,2(3,2MB 


0.274 ± 


014 


-1 


034 







053 





150 


410 


7/464 


jtatistical+Systematic'^ 






















la,l{3,lMB 


0.266 ± 


016 


-1 


116 







081 





154 


448 


4/467 


la,l/3,2MB 


0.273 ± 


016 


-1 


055 







079 





153 


423 


1/466 


la,2P,lMB 


0.276 ± 


016 


-1 


037 







076 





153 


425 


3/466 


2a,l/3,lMB 


0.267 ±0 


016 


-1 


111 







081 





153 


445 


9/466 


la,2l3,2MB 


0.274 ± 


016 


-1 


040 







075 





151 


405 


4/465 


2a,l(3,2MB 


0.272 ± 


016 


-1 


054 







079 





153 


423 


1/465 


2a,2f3,lMB 


0.275 ± 


016 


-1 


041 







076 





152 


421 


1/465 


2a,2p,2MB 


0.274 ± 


016 


-1 


040 







075 





150 


405 


3/464 



'^We exclude the host galaxy systematic term in these fits; see 



5.3 



for details. 



